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Using the improved U(l) and SU(2) slave-boson approaches of the t-J Hamiltonian [Phys. Rev. 
B 64, 052501 (2001)] that we developed recently, we study the doping and temperature dependence 
of superfluid weight and spectral function and discuss our finding of the universal scaling behavior 
of pseudogap. It is shown that at low hole doping concentrations x and at low temperatures T 
there exists a propensity of linear decrease in the superfluid weight n s /m* with temperature T, 
and a tendency of doping independence in the linearly decreasing slopes of -^(x,T) with T in 
qualitative agreement with the experimentally observed relation, that is, -^fir(x,T) = ■^■(x,0) — aT 
with a, a constant. It is also demonstrated that there exists the boomerang behavior, that is, 
both T c and n s /m* increase with hole doping concentration x in the underdoped region, reaches a 
saturation(maximum) at a hole doping above optimal doping and decreases beyond the saturation 
point in the overdoped region in agreement with /xSR measurements. Based on our improved SU(2) 
slave-boson approach to the t-J Hamiltonian, we further investigate the doping and temperature 
dependence of spectral function and discuss our theoretical finding of a scaling behavior of pseu- 
dogap. In addition we discuss the cause of hump and quasi-particle peak in the observed spectral 
functions of high T c cuprates. It is demonstrated that the sharpening of the observed quasi-particle 
peak below T c is attributed to the bose condensation of holon pair in agreement with observations. 
From the predicted ratios of pseudogap Aq to both the superconducting temperature T c and the 
pseudogap temperature T* i.e., fc 2A j! and fc 2/ ^?„ respectively as a function of hole doping concentra- 
tion x, we find that there exists a universal scaling behavior (sample independence) of these ratios 
with doping, by showing a nonlinearly decreasing behavior of the former, that is, £ j ~ x~ a with 

a ~ 2 and a near doping independence of the latter, that is, fc ^?» ~ 4 ~ 6. 



I. INTRODUCTION 



Presently there exist two outstanding physical prob- 
lems related to superfluid weights for high T c cuprates. 
One is the understanding of the so-called 'boomerang' 
(reflex) behavior and the other, the doping independence 
of linear decrease of superfluid weight ^r(x, T) vs. tem- 
perature T. More than a decade ago the transverse field 
muon-spin-relaxation(^-SR) measurements of the mag- 
netic penetration depth A in high T c copper oxide super- 
conductors revealed the 'boomerang' behavior of the su- 
perconducting temperature T c and the superfluid weight 
n s /m*, that is, a linear increase of both T c and n s /m* in 
the underdoped region, a saturation of T c near optimal 
doping and a decrease (reflex) of both T c and n s /m* in 
the overdoped region as the hole doping concentration 
increases^, 0, 0- Earlier Emery and Kivelson^jj sug- 
gested that the linear proportionality of T c and n s /m* 
in the underdoped region can be understood in terms 
of the phase fluctuations of Cooper pairs. It is argued 
that phase fluctuations become important in the under- 
doped region because the superconducting charge carrier 
density increases linearly with hole concentration^. The 
linear increase of the superfluid weight in the underdoped 



region was predicted from both the slave-boson theories 
of the t-J Hamiltonian^, El an d the numerical 

studies of projected trial wave functional]. In the single 
holon slave-boson theories T c is given by the single holon 
bose condensation temperature in the underdoped region 
and the spin gap temperature in the overdoped region. 
Similarly, the Ioffe-Larkin formula^^ shows that the su- 
perfluid weight is mainly determined by the holon contri- 
bution in the underdoped region and by the spinon con- 
tribution in the overdoped region 7J. Thus both T c and 
n s /m* are claimed to increase with hole concentration x 
in the underdoped region and decrease in the overdoped 
region. However the predicted T c increases too stiffly 
with x in the underdoped region because T~ is propor- 
tional to a large hopping energy t ~ 0.5eF[l3} according 
to these single-holon condensation theories. In the non- 
linear sigma model description of the SU(2) slave-boson 
theory [Talll collective modes associated with order pa- 
rameter fluctuations are shown to play important role 
in determining T c and n s /m*. From the projected trial 
wavefunction studies |i"H it was predicted that the super- 
conducting order parameter is arch-shaped as a function 
of hole concentration and that the upper bound of the 
superfluid weight increases linearly with hole concentra- 
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tion. However this study is limited to the T = OK and 
thus the relation between T c and n s /m* has not been 
directly addressed. 

Most recently, magnetic penetration depth(A) 
measurements [T5L Il6l Il7) showed a doping indepen- 
dence in the negative slopes of the superfluid weight 
with increasing temperature showing the relations 
^r(x,T) — ^r(x,Q) — aT in the underdoped region 
where x is the hole concentration and a, a temperature 
independent constant. The observed ratio between 
the superfluid weight and the superconducting tem- 
perature ao = '"t ^')°' > is nearly doping independent 
in the underdoped regionQ. As a consequence the 
normalized superfluid weight displays a universal be- 
havior, following the relation "'/^.^'q) = 1 — P(t^) 
with j3 = -^fiH, Il6l Il7| . This universal behavior has 
been interpreted by the phase fluctuation model (l8|. 
or the quasiparticle model based on the Landau-Fermi 
liquid theory 0, |2(], 0, l^]. However, apart from 
these phenomcnological models the microscopic expla- 
nation of this universal behavior from a microscopic 
theory is unresolved 23|. Earlier other U(l) slave-boson 
theories H S S IH H|| of the t — J Hamiltonian 
predicted that at small hole doping the superfluid 
weight ^p-(x,T) as a function of temperature is strongly 
doping (x 2 ) dependent. This is in direct conflict with 
the /iSR measurements [H E E3- Lee and Wen[j| ^ 
proposed that the SU(2) theory which incorporates the 
low energy phase fluctuations of order parameters may 
resolve this problem. Recently D.-H. Lee[24| presented a 
U(l) slave-boson theory concerned with the low energy 
fluctuations of order parameters and the excitations of 
massless Dirac fermions at the d-wave nodal point. This 
theory, also, showed the x 2 dependence of at finite 
temperature. Most recently Wang et al.[25( also showed 
that the recently proposed d-density wave theory [26) 
fails to predict the doping independence of the linearly 
decreasing slope of superfluid density vs. temperature. 

Angle resolved photoemission spectroscopy(ARPES) 
measurements of high T c cuprates revealed both the tem- 
perature and doping dependence of both spectral peak in- 
tensity and pseudogap at varying pseudogap temperature 
T* and superconducting temperature T c . These mea- 
surements have shown a continuously increasing trend of 
pseudogap (leading edge gap) with decreasing tempera- 
ture even below T c and the appearance of sharp peaks 
below T c [2l|2|. Earlier, Wen and Lee@ reported the 
momentum dependence of the spectral function based 
on their SU(2) slave boson theory involving single-holon 
bose condensation. Chubukov et al. obtained the peak- 
dip-hump feature of the spectral function by using their 
spin fcrmion model 29] . Most recently, Muthukumar et 
al. studied the momentum dependence of the spectral 
function in the resonating- valence-bond state[3(J. How- 
ever there exists lack of comprehensive investigations 
concerning both the temperature and doping dependence 
of spectral functions. 



In this study we plan to pay attention to these prob- 
lems based on our earlier SU(2) slave-boson theory [3l| 
which differs from others in that coupling of the spinon 
pairing order to the holon pairing order A b (holon- 
pair channels) is essential for the formation of Cooper 
pairs as is well depicted in the last term in Eq.J2J), that is, 
.4V :j lAf/fA^^ + c.c.]. Here we would 
like to stress that, obviously, the Cooper pair is a compos- 
ite particle made of the spinon pair and holon pair in the 
slave-boson langua ge. Based on this improved treatment 
over a previous one|32j] we were able to successfully repro- 
duce the arch-shaped superconducting temperature curve 
in the phase diagram of high T c cuprates in agreement 
with observation. Further, the application of this theory 
resulted in the peak-dip-hump structure of optical con- 
ductivity again in agreement with observations j3^. Us- 
ing this theory we present a five-fold concerted study on 

1. the temperature dependence of the superfluid weight, 

2. the boomerang behavior of the superfluid weight and 
the superconducting temperature, 3. the origin of the 
hump and the sharp quasiparticle peak in observed spec- 
tral functions (ARPES)^ 4. the doping and tem- 
perature dependence of the spectral functions, and 5. the 
universal (that is, sample independent) scaling behavior 
in the ratio of the spin gap in the spectral function to the 
superconducting temperature and the pseudogap tem- 
perature respectively with hole doping concentration, all 
of which are found to originate from the coupling effect 
of the spin (spinon pairing order) to the charge (holon 
pairing order) degrees of freedom. 



II. TEMPERATURE DEPENDENCE OF THE 
HOLON AND SPINON PAIRING ORDER 
PARAMETERS 

Earlier, we presented an improved U(l) and SU(2) 
slave-boson approach of the t- J Hamiltonian which differs 
from other slave-boson approaches to the t- J Hamiltonian 
in that coupling between the charge and spin degrees of 
freedom is manifested in the slave-boson representation 
of the Heisenberg term of J X]<i j> (^i ■ Sj — \niiij). It 
is noted that the charge degree of freedom is well exhib- 
ited in the intersite charge interaction term of ^n.n-, [3l|. 
Owing to the coupling of the holon pairs to the spinon 
singlet pairs as mentioned in the previous section, the 
Cooper pairs are formed as composite particles of the 
spinon pairs and holon pairs. Holons form holon pairs via 
coupling to spinon pairs and the holon-pair bosons un- 
dergo bose condensation 35] at and below T c . This is the 
major difference of the present holon-pair boson theory as 
compared to other earlier theories which pay attention to 
singe-holon bose condensation^ 0, B B liil I2i4 HH • 0m 
theory was shown to reproduce the arch shaped super- 
conducting (holon-pair bose condensation) temperature 
T c as a function of doped hole concentration 2; [3]]], in 
agreement with the experimentally observed phase dia- 
gram. It is, thus, of great interest to employ this theory 
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to further check its applicability to other important ar- 
eas of physics; the doping dependence of the superfluid 
weight and the spectral function. Here we focus our at- 
tention to the doping and temperature dependence of the 
order parameters. 

Since the predicted results of various physical proper- 
ties are qualitatively indistinguishable between the U(l) 
and SU(2) slave-boson theories, below we present only 
the SU(2) theory for generality. In the SU(2) slave- 
boson representation^, the electron operator is given 

/i 



and the spinon action, 



by c a = -^h^ipa with a 



1,2, where ipi 



h 



blets of spinon and holon annihilation operators in the 
SU(2) theory. The SU(2) slave-boson representation of 
the t-J Hamiltonian shows |3l| 

<i,j> 

+(/ 2 i/y - /li/jyX&y&M + b\ t b 2j ) + h.c. 

~\ E ^ + h\h l )(l + h]h J ) x 

<i,j> 

iflifli flift^ihihi - fa hi) - i>Y. h >- 

i 



and h 



are respectively the dou- 



(1) 



where ^W'^M 3 ) are rea j Lagrangian multipliers to 
enforce the local single occupancy constraint in the SU(2) 
slave-boson representation It is noted that the in- 
clusion of the charge-charge interaction —jriiTij in the 
Heisenberg Hamiltonian results in the coupling between 
the spin (spinon) and charge (holon) degrees of freedom. 
It is easy to show that under inverse SU(2) transforma- 
tion of the above expression we can readily recover the 
original form of the U(l) slave-boson representation (see 
|3ll | for details). 

After relevant Hubbard-Stratonovich 

transformations |3lj . the holon action is obtained to 
be 

S b ( A, x , A y , A B , h) = f dr \ V {v u r) (d T - 
Jo i 

4 E K-i 2 (Ei A Ui 2 

<i,j> a.J3 
<i,3> 

J2 (|A£ffct (r . )T)A ^t (r . )T)) 



S f (x,A f ,i>) = f dr\y^^{v h r)d T ^{v U T) 
Jo t 

I J(l-^ 2 ) £ (JA'' 2 ' 1 '--' a 



<i,j> 



Here various symbol definitions are as follows. ft,(r,-,r) = 
h(ri,r) \ ig ^ doublet of holon field; and 



b2{ri,T) 



/i(""i,r) 



the SU(2) doublet of spinon 



field. A is the electromagnetic (EM) field. U\ a 



X*. -A f - 
-A* 



and U, 



f 

id 



-2A{* 



- 2 4 



are 



-y -Xij ) \ ~ liX h ~Xv ) 

the order parameter matrices involved with hopping(xy) 



' 'Miii^ + 0/(r;)), with r] 



and spinon pairing(A^ ). Here Xi,i+i — TJ e '"' ' ' ' ' 



flt fo)) and A{ l+l = tie- 
yj\x\ 2 + \Af\ 2 and 6»j° = ± tan" 1 ^-(thc sign +(-) is for 
the ij link parallel to x (y)). a<;(rj) represents the phase 
fluctuations of the hopping order parameter, /3/(rj), the 
phase fluctuations of the spinon pairing order parame- 
ter, and #/(i"i), the relative phase fluctuations between 
the amplitudes of hopping and spinon pairing order pa- 

{ A b A b - \ 
rameters. A^j = « aI 3 ' 12 I is the matrix of the 

V V A ' ,:2I A z,;22 / 

holon pairing order parameter. 

After integration over the holon and spinon fields, we 
obtain the total free energy, 



F(A) = ~1n J D x DA f DA B 

-p(F b (A,xA f A B )+F f (xA f )) 



(4) 



+ c.c 



J 
2 



+ c.c. 



<h3> 



where F b (A, x, A? ,A B ) = - \ In J Dhe- s "( A ^ Af ' AB ^ 
is the holon free energy and Ff( x ,A f ) 
— 4 In J Dijje~ sl ( x > Al the spinon free energy. 

The saddle point equations are solved for the spinon 
pairing order parameter A f and holon pairing order pa- 
rameter Af, including the hopping order parameter as a 
function of doping and temperature. In Fig. 1 we dis- 

, r)play the temperature dependence of the order parame- 
ters A f and A 6 at x = 0.1 for J/t = 0.2, 0.3 and 0.4. 
The spin gap temperature T* and the superconducting 
temperature T c are identified as the temperatures below 
which the spinon pairing and holon pairing gaps respec- 
tively begins to open up. Both the spinon pairing and 
the holon pairing order parameters are shown to change 
continuously at the critical temperatures, T* and T c . In 

(2)Fig. 2 we show the doping dependence of T c and T* 
for J/t = 0.2, 0.3 and 0.4. The spin gap temperature 
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T* is predicted to increase with increasing antiferromag- 
netic coupling J as expected. It is noted that T c also 
increases with increasing J. This is because the Cooper 
pairs as composites of the holon pair and the spinon pair 
are bose condensed as a result of coupling of the holon 
pair to the spinon pair whose order parameter strength 
increases with J, as can be understood from the last term 
in Eq. @. This is in sharp contrast to the single holon 
slave-boson theories00,|3,l3 where T c scales only with 
t, but not with J. In Fig. 3 the doping dependence of 
the spin gap 

A = J(l — x) 2 Af(cosk x — cosky) (5) 

at k — (tt, 0) and T = is shown for the choice of J/t — 



0.2, 0.3 and 0.4. The spin gap monotonically decreases 
with increasing hole concentration. This is because the 
spin (spinon pair) singlet bonds are more readily broken 
by mobile holes with increasing vacant sites. 



III. DOPING AND TEMPERATURE 
DEPENDENCE OF SUPERFLUID WEIGHT 

To compute the EM current response function, we first 
obtain the total free energy by integrating out the three 
phase fields a, f} and 9 above. We obtain a formula for 
the EM current response function, up to the second order, 



IL m („,q) = ILT' V,q) 

- E E n^'i^^n^^j + n/i^C^n^Kq), (6) 

a 1 ,a 2 =e,a,/3 l' ,m'=x,y 



where [U b (co, q) + IP>, q)] J /[q . 
trix element of 



is the inverse ma- 



IT 



b(a\a 2 ) 



and II 



nb(a ,a ) I \ 
f(a\a 2 ) 



I m 



>,q) 



Here 



represent the isospin current re- 



sponse functions of holons and spinons respectively to 
gauge fields a 1 , a 2 = 6, a, or /3. Similarly, n^'^ (uj, q) 
is the EM current response function of holon to the EM 

field. n?^f ' A ' is the isospin current response function 
of holon to both gauge fields a and A. It is noted that 
the SU(2) isospin current of holon(spinon) is defined as 

6(/) - -/3ff^-, and the 'EM' current, j t = -fijj-. 



.) 



<\.i 



In the SU(2) slave-boson theory, the response function 



IT 



fe(A,a 1 ) 
I' 



of the holon isospin current to the EM field van- 



ishes owing to the contribution of the 62-boson in the 
static and long- wavelength limit 9] . Therefore the super- 
fluid weight of the total system is given only by the holon 
current response function, 



n s 
m * 



-— lim I!,, 

e 2 q-+o " 



b(A,A) 



(w = 0,q) 



(7) 



Here 11^^ (lu, q) is computed from the use of the usual 
linear response theory for the holon action Eq. (J2J) 36J . 

In Fig. 4(a), with the use of our slave-boson theory[3lJ 
we show the temperature dependence of superfluid weight 
for a wide range of doping concentrations covering both 
underdoped and overdoped regions. The slope of each 
curve represents the variation of superfluid weight with 
temperature which is only in qualitative (bu t not quan- 
titative) agreement with observations [l6l Il7j in that 
the superfluid weight -^(£,0) at T = OK is overesti- 
mated and further -^(x, T) reaches zero rather abruptly 



compared to the observation which showed a slower 
drop0|. However a close look at Fig. 4 for comparison 
with experiment 0] reveals a linear slope behavior only 
in the lower temperature region and a quick non-linear 
drop of superfluid weight at the higher temperature re- 
gion. Such a slope variation is in good agreement with 
our results although our theory predicted a faster drop 
in -^lr(T). The continuous drop in the predicted super- 
fluid weight with temperature indicates a second order 
phase transition. We believe that improvement can be 
achieved if we can introduce the t term in our t- J Hamil- 
tonian which is known to cause a good agreement with 
the observed dispersion energy by showing realistic cold 
and hot spots. The contribution of the hot spot in the 
Brillouin zone becomes increasingly important with in- 
creasing temperature for classical phase fluctuations to 
occur at T c . 

As temperature T increases, the predicted superfluid 
weight shows a tendency of linear decrease in T with 
nearly identical slopes particularly at low temperatures 
in the underdoped region. To rigorously check such a 
propensity of doping independence in the slope of 
vs. T, we computed differences in the superfluid weight 
at T = and T ^ 0, that is, %fir{x,T) - ^t(x,0) at 
each x. The results are shown in Fig. 4 (b). A de- 
crease of the superfluid weight with temperature at all 
hole concentrations is predicted although the tempera- 
ture dependence of the superfluid weight does not show 
a clear linearity at all temperatures as shown in Fig. 
4 (b). This figure displays that in agreement with ob- 
servation there is a tendency of linear decrease with a 
nearly identical slope particularly in the underdoped re- 
gion only at low temperatures. Although the temper- 
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ature dependence of ^(x,T) does not show a global 
linearity by yielding a good fit to the empirical relation 
of %t(x,T) =p(x,0) - a(x)T - a(x)T 2 - 0(T 3 ), the 
linear slope a is predicted to be doping independent only 
in the underdoped region and at low temperatures as is 
clearly manifested in Fig. 4 (b). As shown in the figure 
it is found that the variation of the slope, that is, the 
superfluid weight n s /m* vs. T/t at low temperatures 
does not appreciably change in the underdoped region 
(x = 0.04, 0.07 and 0.1). On the other hand, in the over- 
doped region (x = 0.16, 0.18 and 0.2) the slope of the 
superfluid weight with temperature is no longer doping- 
independent quite unlike the underdoped case, again, in 
agreement with observation. The predicted optimal dop- 
ing is x — 0.13. We find that there exists crossings 
of superfluid weights between the underdoped and over- 
doped cases in the plane of n s /m* vs. T/t as shown 
in Fig. 4(a). This prediction is also in agreement with 
measurements 1 1 5l I16L Il7| . 

In Fig. 4(c) we display the the normalized super- 
fluid weight ( " s / m }, x ''^} ) as a function of the normalized 
temperature (T/T c ) for various hole concentrations. As 
shown in the figure, the normalized superfluid weight be- 
gins to decrease linearly at low temperature and drops 
more rapidly at higher temperature, by showing a con- 
vex shape in agreement with experiments Il6t Il7| . 
The predicted decrease of superfluid weight with tem- 
perature occurs as a consequence of diminishing spin 
singlet pairing order with increasing temperature. This 
is readily understood from the holon pairing term (the 
last term in Eq. (2)) which reveals the coupling of the 
spin (spinon) singlet pairing order with the holon pairing. 
Consequently the decrease of the superfluid weight with 
increasing temperature is caused by breaking the bond 
of spinon singlet pairs coupled to the holon pairs, whose 
spin bond strength decreases as temperature increases. 
Although the superfluid weight shows a linear decrease 
in agreement with observation, the coefficient (3 in the 
relation of ""^/".^'T? = 1 — d(S-) is scattered around 

n s /m* (£,0) "^T c ' 

0.055 and thus fails to quantitatively agree with the em- 
pirical value of /3 ~ 0.5. However the important finding 
here is the (3 is independent of doping in the underdoped 
region 17]. The quantitative discrepancy indicates that 
the inclusion of the next nearest hopping (t ) term is nec- 
essary. This is because the Fermi surface segment grows 
from the nodal point initially at half-filling and closes 
near (tt, 0) as hole doping reaches optimal doping and 
beyond as is well-known from the an gle resolved photoc- 
mission spectroscopy measurements [37j . In the future 
the inclusion of the t term will have to be made to take 
a good account of classical phase fluctuations in the re- 
gion of the hot spot at high temperatures at and near 
T c 

According to the present theory the superconductiv- 
ity arises by the condensation of the Cooper pairs of d- 
wave symmetry as composites of the holon pairs of s-wave 
symmetry and the spinon pairs of d-wave symmetry, but 



not by the single-holon bose condensation unlike other 
theories. To put it otherwise the Cooper pairs can be 
regarded as composite objects resulting from such cou- 
pling of the spin (spinon) and charge (holon) degrees of 
freedom. 

To clearly show the effect of coupling between the 
spin(spinon) and charge(holon) degrees of freedom, we 
calculated the superfluid weight by varying the antifer- 
romagnetic spin-spin coupling, that is, the Heisenberg 
coupling J. In Fig. 5, we show the J dependence of su- 
perfluid weight at x = 0.1 for three different choices of 
J/t = 0.2,0.3 and 0.4. As J increases, a decreasing ten- 
dency in the slope of the superfluid weight vs. temper- 
ature is predicted. For large antiferromagnetic coupling 
J and thus the strong spinon pairing bond (strength), it 
is more difficult to thermally break the holon-pairs be- 
cause of the increased holon-spinon coupling. It is again 
reminded that the Cooper pair is obviously a composite 
of the spinon pair and the holon pairs which are coupled 
to each other as is shown in the last term of Eq. (01 • This 
results in the more slowly decreasing slope of the super- 
fluid weight with increasing temperature for large J, as 
is shown in the figure. 

In Fig. 6, we show the predicted boomerang behavior 
in the plane of T c vs. -^(x, T — > 0), by showing a linear 
relationship between the two at low hole concentrations 
in the underdoped region and the 'reflex' behavior in the 
overdoped region based on the U(l) theory. By the re- 
flex behavior we mean the decrease of both T c and 
as the hole doping concentration increases in the over- 
doped region. This predicted boomerang behavior is con- 
sistent with the measurements pi 0, Q of the muon-spin- 
relaxation rates a (a oc n s /m*). Here we focus our atten- 
tion to its cause. The pseudogap (spin gap) temperature 
decreases in the overdoped region in agreement with ob- 
servation. The predicted pseudogap (spin gap) caused 
by the appearance of the spinon pairing order A-^ de- 
creases in the overdoped region and this, in turn, causes 
a decrease in holon pairing A b owing to its coupling to 
the spinon pairing order whose strength decreases with 
increasing hole concentration. As a consequence, has 
to show the boomerang behavior in the plane of vs. 
T c . The boomerang behavior of the SU(2) theory is found 
to be qualitatively the same as that of the U(l) theory. 



IV. UNIVERSAL SCALING BEHAVIOR OF 
PSEUDOGAP BASED ON THE HOLON-PAIR 
BOSE CONDENSATION THEORY 

The single particle(electron) propagator of interest is 
given by a convolution integral of spinon and holon prop- 
agators in the momentum space jj-j, 

G a0 (k,u;) = % - y ^^[^G^k+k'.w + w') x 

l,m 

G^(kV)]- (8) 
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The mean field Green's functions are G^; m (k, w) = 

-ijdt^e™*-*-* < T[^(x,f)^ m (0,0)] > and 

G? m (k, w ) = -i/dtE.e*-*"*-* < r[6 I (x J *)6j n (0,Q)] > 
respectively (The symbol < > here refers to the finite 
temperature ensemble average of an observable quantity 
O, < O >= itr(e-^O)). 

The one electron removal spectral function, 7(k, u>) is 
obtained from [ssj . 

I(k,w) = -ilmG(k,w + iO+)/(a;), (9) 

7T 

where /(w) is the Fermi distribution function. The 
Heisenberg coupling constant of J = 0.2 t and the hop- 
ping strength of t = 0.44 eV"12| are chosen in the 
present calculations. Using the SU(2) theory, the pre- 
dicted values of optimal hole doping x Q , pseudogap tem- 
perature T* and bose condensation temperature T c are 
x Q = 0.13, T* = 0.029t(U8K) and T c = 0.021t (107.27T) 
respectively. For the evaluation of the spectral functions 
7(k, lo), the line width of the Lorentzian function is set 
to be e = 0.01i(4.4mey) with the choice of t = 0.44eV. 

It is known from the ARPES measurements of 
LSCOlU and BSCCO[I3 that a g° od Fermi nesting 
appears in a region near the M point (k = (ir,0)) at 
which maximal antiferromagnetic correlations between 
electrons are realized. The predicted binding energy 
monotonically decreases with the increase of momentum 
k along T - M, in agreement with the ARPES |4~l|. 

In Fig. 7 we display the computed spectral functions 
7(k, u>) at k = (w, 0) both below and above T c . The 
computed gap size at k = (w, 0) is predicted to coin- 
cide with the spin gap size A in Eq.JSJ). The peak 
position occurs at smaller binding energies at temper- 
atures (0.022t(112.3iT)) above T c . As is shown in Fig. 
7(a), a spin gap begins to open at a higher critical 
temperature T* = 0.029t(147.97C) above T c . The pre- 
dicted peak position shift, that is, the shift of lead- 
ing edge gap is much slower below T c . This trend is 
in complete agreement with observations 03- 
cently ARPES measurements^ for Bi 2 Sr2CaCu20s 
showed that there exists a quasiparticle peak above T c 
and the peak intensity decreases with increasing temper- 
ature, while the enhancement of peak intensity occurs 
markedly at temperatures below T c . Our computed re- 
sults revealed the presence of the quasiparticle peaks in 
both below and above T c as shown in Fig. 7 (a). A 
recent ARPES study of Feng et al. Q claimed that the 
enhancement of the quasiparticle peak below T c is caused 
by bose condensation. Our theory proves that such en- 
hancement of the quasiparticle peak below T c is caused 
by the bose condensation. For this verification, in our 
calculation we first removed completely the contribution 
of the holon-pair channels (which comes from momenta 
k = (0, 0) and k = (it, it)). The sharp peak is now seen 
to disappear and only the hump structure is predicted 
at all temperatures at and below T c while the pseudo- 
gap size remains unchanged. This clearly indicates that 



the pseudogap is caused by the formation of spin sin- 
glet pairs and that the presence of holon-pair bosons is 
essential for yielding the observed enhancement of quasi- 
particle peaks at the M point below T c . In short, the 
quasiparticle peak is caused by bose condensation below 
T c . We may now argue that the appearance of the hump 
in the absence of bose condensation is seen to be caused 
by the antiferromagnetic spin fluctuations of containing 
the shortest possible correlation length, that is, the spin 
singlet pair excitations. Both the SU(2) and U(l) slave- 
boson theories were found to predict the enhancement 
of quasiparticle peaks at k = (it, 0) and k = (tt/2,tt/2) 
below T c . As is shown in Fig. 7, the predicted leading 
edge gap at k = (tt, 0) showed a continuous increase as 
temperature decreases from a pseudogap temperature T* 
to temperatures below the superconducting temperature 
T c . This indicates that the observed leading edge gap 
with the appearance of sharp quasiparticle peak below 
the superconducting transition temperature is originated 
from the persistence of the spin gap with the emergence of 
bose condensation (holon pair bose condensation). This 
implies that the observed quasiparticle peak in the su- 
perconducting state is attributed to the presence of the 
coupling of the spin (spinon pairing order A-^) to the 
charge (holon pairing order A b ) degrees of freedom as 
shown in the last term of Eq.|J2J. 

In Fig. 8 the doping dependence of spectral functions 
is displayed for underdoped, optimally doped and over- 
doped cases at a low temperature, T — 0.004i (T = 
20 AK) below T c , near the M point (k = (0.8tt,0)) 
for comparison with observation. In agreement with 
the ARPES [Hi the predicted spectral weight of the 
sharp quasiparticle peaks is seen to increase as the hole 
concentration increases upto a tested value in the over- 
doped region. The leading edge gap is shown to de- 
crease with increasing hole concentration in agreement 
with observations |39j . In our calculations the dip in the 
peak-dip-hump structure is not predicted. It is claimed 
that the dip is caused by electron-phonon coupling 46] , 
which is not considered in our present treatment of the 
t-J Hamiltonian. 

Figs. 9. (a) displays as a function of doping rate the 
ratios of the pseudogap 2A (fc = (tt, 0)) at OK to the su- 
perconducting temperature T c and the spin gap temper- 
ature T*, that is, ^=£- and fc 2Z ^, respectively. It is pos- 
sible that each high T c cuprate(e.g., LSCO, YBCO and 
BSCCO) may have an effectively different J value caus- 
ing variation in T c . It is quite encouraging to find that 
in agreement with observations [47] . the ratio 2A °^"^ 7r '°^ 
decreases rapidly with hole concentration x, showing a J 
independent and thus sample- independent universal scal- 
ing behavior with x. On the other hand, 2Ao< k^T*'°^ 
shows a nearly doping independence of ranging its value 
between 4 and 6, as shown in Fig. 9.(b). The doping 
independence of 2A °0*-(Tr ,o)) - g ^ ue tQ ^ e f ac t that both 
Ao and T* simultaneously decreases with i in a con- 
certed manner as is shown in Figs. 2 and 3. That is, the 
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ratio should remain unchanged by revealing 

that both Ao and T* decrease at an equal rate with in- 
creasing hole concentration x. On the other hand, the 
non-monotonously decreasing behavior of 2 Ao(fc— j>,0)) is 
attributed to the arch shape of T c line as a function of 
hole concentration^. The arch shape bose condensa- 
tion temperature is due to the coupling of the holon pairs 
and spinon pairs. In the underdoped region, spinon pair- 
ing order A* is large and T c increases with increasing x. 
On the other hand, in the overdoped region continu- 
ously decreases with increasing x. In this region T c also 
decreases with increasing i as a result of the diminish- 
ing nature of A f . It is found that and do 
not appreciably change with the variation of Heisenberg 
coupling J. This is due to the nature that Ao, T* and 
T c are found to be proportional to J, keeping these ra- 
tios unchanged Thus we found that independent of 
Heisenberg exchange (antiferromagnetic) coupling there 
exists a universal scaling behavior of the ratios of Ao /T c 
and A /T* respectively. 

V. SUMMARY 

We stress that the Cooper pair is simply a compos- 
ite of spinon pair and holon pair as a consequence of the 
mutual coupling in the slave-boson language. The spinon 
pair (spin singlet pair) formation occurs below the spin 
gap (pseudogap) temperature. Thus it is needless to say 
that bose condensation or superfluidity should occurs be- 
low the spin gap temperature. From the present study we 
found the correct doping and temperature dependences of 
both the superfluid weight and spectral function and the 
universal scaling behavior of pseudogap in agreements 
with observations. Regarding the doping and tempera- 
ture dependences of the superfluid weight, we observed 
the nearly doping independence in the negative slope of 
superfluid weight vs. temperature in the underdoped re- 
gion, thus correctly reproducing the experimentally ob- 
served relation, %t(x, T) — ^fir(x, 0) — aT with a, inde- 
pendent of hole concentration. Our results which showed 



a linear decrease of the superfluid weight are qualitative, 
but not numerically accurate as pointed out earlier. On 
the other hand, in the overdoped region the predicted 
slope of the superfluid weight varies with both hole con- 
centration and temperature in agreement with observa- 
tion. The boomerang behavior in the locus of both 7% 
and T c with the variation of hole concentration in the 
plane of ^(x : T — > 0) vs. T c is in agreement with the 
/i-SR experiments. The decreasing (boomerang) behav- 
ior of j^r and T c in the overdoped region is attributed to 
the weakened spinon pairing order coupled to the holon 
pairing order. From the present study we find additional 
salient features in agreements with observations. They 
are 1. the sharp quasiparticle peak observed below T c 
is attributed to the bose condensation, and the hump is 
caused by the antiferromagnetic spin fluctuations of con- 
taining the shortest possible correlation length namely 
the spinon (spin) singlet pair excitations, 2. the spectral 
intensity of quasiparticle peak shows an increasing trend 
with hole concentrations in the underdoped region and 
the pseudogap (spin gap) size at k = (tt, 0) continuously 
increases as temperature decreases from the pseudogap 
temperature T* to temperatures at and below the super- 
conducting temperature T c , 3. a decreasing trend of spec- 
tral peak intensity is predicted as temperature increases, 
and, 4. there exists a universal (i.e., sample indepen- 
dent) scaling behavior of both and fc 2 ^?, with x. 
Finally coupling between the spin (spinon pairing order 
A^) and the charge (holon pairing order A b ) is impor- 
tant to explain the observed phase diagram, the doping 
and temperature dependence of the superfluid weight, its 
boomerang behavior, the doping and temperature depen- 
dence of the spectral function and the universal scaling 
behavior of the pseudogap. We stress again the Cooper 
pair is a composite of the spinon pair and the holon pair 
to yield an object of charge +2e (or — 2e) and spin 0. 
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FIGURE CAPTIONS 

Fig. 1 Temperature dependence of the spinon pairing or- 
der parameter (A/) and the holon pairing order 
parameter (A&) at x — 0.1 for J/t = 0.2, 0.3 and 
0.4. 

Fig. 2 Doping dependence of the superconducting tem- 
perature T c and the spin gap temperature T* for 
J/t = 0.2, 0.3 and 0.4. 

Fig. 3 Doping dependence of the spin gap Ao(fc = (ir, 0)) 
at T = for J/t = 0.2, 0.3 and 0.4. 

Fig. 4 (a) Temperature dependence of the superfiuid 
weight Zfc for underdoped(a: = 0.04,0.07,0.1), 
optimal doping(a; = 0.13) and overdoped(x = 
0.16,0.18,0.2) rates with J/t = 0.2. (b) The dif- 
ference of superfiuid weights between T ^ OK and 
T = OK, that is, %k(T)-2fir(T = 0). Each line rep- 
resents a fitting to the computed superfiuid weight, 
(c) The normalized superfiuid weight „ s /m* (z'o) 
as a function of scaled temperature T/T c for same 
doping concentrations and J value as (a). 

Fig. 5 The difference of superfiuid weights between T 7^ 
OK and T = OK, £fc(T) - ^.(T = 0) for a un- 
derdoped hole concentration x = 0.1 with J/t — 
0.2, 0.3 and 0.4. Each line represents a fitting to 
the computed superfiuid weight. 

Fig. 6 The superfiuid weight j^(x,T) at T = O.OOlt 
(equivalent to T ~ 5K with the use of t = 
0.44eV[T3') vs. superconducting temperature T c 
with the choice of J/t — 0.2 based on the U(l) 



theory. The open box represents hole concentration 
starting from x = 0.01 to x = 0.22 and the ar- 
row denotes the direction of increasing doping rate. 
The predicted optimal doping rate is x = 0.07 with 
the U(l) theory. The SU(2) theory predicts quali- 
tatively the same boomerang behavior. 

Fig. 7 Temperature dependence of the spectral functions 
/(k, uS) for k = (tt, 0) at the predicted optimal 
doping of x = 0.13 (T c = 0.021t(107.2Jr), T* = 
0.029t(147.9iC)) with J = 0.2t for k = (tt,0). (a) 
/(k, ui) with all momenta of holon included and (b) 

7(k, w) with the exclusion of only k = (0,0) and 

(7T,7r). 

Fig. 8 Doping dependence of spectral functions 7(k, u>) 
near the M point (k = (0.87T, 0)) at temperature 
0.00At(20AK) below T c for underdoped(z = 0.04 
and x — 0.085), optimal doping(a; = 0.13) and 
overdoped(a; = 0.2) rates. 

Fig. 9 The ratios of the spinon pairing gap at T = OK to 
(a) the superconducting temperature, 2A °^ ; ^ 7r ' ^ 
and (b) the spin gap temperature, 2A o(^-An,o)) ag 
a function of j/- where x = 0.13 is the predicted 
optimal doping concentration for J/t — 0.2 and 0.3, 
and x = 0.14 for J/t = 0.4. The experimentally 
obtained universal ratios for Bi2Sr2CaCu20s+x, 
La2- x Sr x Cu04, and ThBa2CuO^ are denoted as 
solid circles |47|| . 
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